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In practically all turbulent flows, turbulent energy decay is present and competes with numerous 
other phenomena. In Kolmogorov's theory, decay proceeds by transfer from large energy-containing 
scales towards small viscous scales through the "inertial cascade." Yet, this description cannot 
predict an actual decay rate, even in the simplest case of homogeneous isotropic turbulence (HIT). 
As empirically observed over 50 years, the steepness of the "infrared" spectrum — at scales larger than 
energy-containing eddies — determines decay, but theoretical understanding is still missing. Here, 
HIT decay laws are derived from angular momentum invariance at large scales — an approach first 
mentioned by Landau in 1944, but unduly dismissed later. This invariance also predicts the so-far 
unexplored turbulence decay in layer, tube, and spot zones in a fluid at rest. Beyond expanded 
and simplified theoretical descriptions, these findings suggest new practical modeling strategies for 
turbulent dissipation, often deficient in applied simulations. 



Background. Incompressible homogeneous isotropic 
turbulence (HIT) at high Reynolds number evolves un- 
der the sole influence of its own dissipation: no external 
source of a characteristic quantity exists (flow size or ve- 
locity, energy production, etc.). After some initial tran- 
sients, it is thus expected to reach a self-similar behavior 
whereby, with the proper origin on time t, 
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k being the (per mass) mean turbulent kinetic energy, 
and A the integral length scale — the typical size of the 
large energy-containing turbulent eddies. As first noted 
by Kolmogorov 1 , self-similarity of decay reflects the ex- 
istence of an invariant in the turbulent field, which can 
be written as 



I = k\ n 



and thus n = 
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according to Eqs. (1). Experiments (wind tunnels, water 
channels, cryogenic helium, etc.) provide values with a 
clear clustering around n ~ 1.2 2,3 corresponding to m w 
3. 

This simple dimensional analysis of mean values is 
not sufficient to characterize neither the underlying ran- 
dom velocity field u(r, £), nor the fluid dynamics process 
which produces an invariant starting from a given ini- 
tial condition u(r,0). The next descriptive order is pro- 
vided by the velocity correlation tensor Ui(0,t)uj(r,t) 
(the overbar stands for ensemble averaging). With the 
more common representation in Fourier space, it can be 
reduced to the spectral density of turbulent kinetic en- 
ergy E(k) = J [J I u t (0) Ul (r) e iK - r d 3 r] K 2 d 2 uj K thanks 
to the conditions of average homogeneity and isotropy, 
and thus k = J E(k)c\k. In fully developed turbulence 
E(k) peaks around k\ sw ir/X (see Fig. 1) and follows Kol- 
mogorov's scaling E(k) ~ k~ 5 / 3 in the inertial range — 
from k\ up to Kolmogorov's dissipation wave number. 

The infrared range of E(k) — below n\ — was very early 
found to be related to the invariant I 4 ' 5 ' 6 ' 7 , but it is diffi- 
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Figure 1: Schematic spectral density profiles of energy for 
HIT (in logarithmic coordinates scaled by the time-dependent 
energy k and the integral length scale A). Depending on its 
logarithmic slope /i, the infrared range (i) transitions quickly 
towards \i = 4 for /i < 4, (ii) slowly evolves towards /i ~ 3 
for 3 < fi < 4, or (iii) is invariant for /i < 3, yielding self- 
similar decay of HIT. The inertial range decays but preserves 
its Kolmogorov —5/3 slope. 



cult to produce and control in experiments. Only recently 
could its influence be thoroughly observed by numerical 
means 8 ' 9 , showing three different regimes depending on 
the steepness of its self-similar profile (see Fig. 1): 

(i) For li > 4, a quick transition towards ix = 4 takes 
place, typically in a few turnover times X/y/k. This well- 
documented process 5 designated as "backscattering" is 
due to the coupling of wave numbers above k\ into elon- 
gated triads through the non-linear term of the Navier- 
Stokes equation. 

(ii) For 3 < fi < 4 — including li — 4 as resulting from 
point (i) above, — a slow, non-self-similar asymptotic evo- 
lution towards li ~ 3 takes place. This produces two 
close but different regimes within the infrared range (see 
Fig. 1), but an almost self-similar evolution of k(t) can 
eventually be observed for late-enough times. 

(iii) For li < 3 — including li ~ 3 as resulting from point 
(ii) above at very long times, — a full self-similar decay 
of the spectrum is obtained, i.e. E{n) is invariant with 
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respect to the reduced units k, A and their combinations. 
Moreover, the infrared profile is fully invariant in or- 
dinary units — this is the so-called "permanence of large 
eddies" — and the overall invariant of turbulent decay is 
thus given by m = // + 1 . The initial preparation of the 
infrared spectrum of HIT thus defines the invariant. As 
explained by Saffman 6 , m = 3 or fj, = 2 prevails in most 
experiments if they all tend to produce randomly dis- 
tributed initial impulses on the energy containing struc- 
tures in the field. 

In summary, it is well understood how an assumed 
permanence of large eddies produces an invariant ' 6 ' 7 — 
ultimately controlling n, — how an initial k 2 infrared pro- 
file is generated 6 , and how an initial k m profile with // > 4 
is quickly reduced to fi = 4 5 . However, a formal proof of 
the permanence of large eddies for [i < 3 does not seem 
to have been produced so far 5 ' 8 ' 9 . The present article is 
focused on this critical aspect and its extension to situa- 
tions other than HIT. 

Incidentally, one should notice that in the present work 
the far infrared and dissipative ranges — at K typically be- 
low or above the ranges represented in Fig. 1 — are not 
constrained to evolve in a self-similar manner, because 
they marginally affect the overall behavior of energy con- 
taining scales. In contrast, the infrared contribution to 
k is by no means negligible: not only does it control the 
decay rate, but it also represents about |/(^ + 1) of the 
energy contained in the inertial range, a substantial 22% 
for \i = 2. 

Landau's invariants. The first and most direct rela- 
tionship between the infrared spectrum and a decay in- 
variant was postulated by Loitsyanskii 4 in direct space, 
using the normalized longitudinal space-correlation func- 
tion of velocity f(s, t) = u z (Q, t) u z (sz, t) / u 2 (0, t). f(s) 
is related to E(n) by Fourier transform, and thus the 
inertial and infrared (or large scales) ranges correspond 
respectively to f(s) w 1 — (s/A) 2 / 3 for s < A, and 



f(s) oc (s/A)" m for s > A. Thus f(s) 
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scales for Saffman's E(k) ~ n 2 infrared spectrum. 

The evolution of f(s,t) is given by the Karman- 
Howarth equation which is deduced from the Navier- 
Stokes equation and involves a normalized two-point 
triple-correlation function K(s). In space integrated form 
it yields the invariance of Loitsyanskii's integral 4 

POO 

hoi.=k f(s,t)s 4 ds, (3a) 
Jo 

cxfc(i)A 5 (0 if f(s,t)=f(s/\(t)), (3b) 

provided that /(s) and K(s) decrease faster than 
s~ 5 . When it was found 5 ' 6 that f(s) and K(s) 
actually decrease slower than s~ 5 in all prac- 
tical situations — as mentioned in the introductory 
paragraphs — Loitsyanskii's approach was simply dis- 
missed. 

Yet, when /Loi. was still considered invariant, 
Landau suggested a particularly appealing physical 



interpretation ' . From dimensional and scaling argu- 
ments, he remarked that Loitsyanskii's integral was sim- 
ilar to the per- volume angular-momentum variance of an 
arbitrarily large spherical volume Q(R), which he esti- 
mated to grow as R 3 : Jlol oc lim HpHp/R 3 (no sum 
on z) as R — > oo, where H = J„ r X ttd 3 r. As he also 
estimated the average torque Tp(R) on i?(i?) to grow 
as R 2 only, the invariance of Jl<h. was recovered. How- 
ever, Landau did not develop his argument in full, and, 
althoug h it appears convincing, it could as well lead to 
HpHp ~ i? 5 instead of R 3 by a similar scaling analysis: 
each point yields a typical kR 2 contribution to Hp Hp 
and the superposition of uncorrelated contributions in 
the double integral yields a further R 3 . 

Landau's interpretation actually reduces to solving the 
stochastic differential equation dfH z 2 (R) = Tp(R) — 
similar to Langevin's equation for Brownian motion — 
and seek invariance conditions. The explicit calculation 
(see Supplementary Methods, §§1-3) eventually yields 
neither of the expected results as 

d t HpHp ~ (A/Vfc) TfTf fori?^oo, 
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where f a are the five independent normalized two-points, 
quadruple space-correlation functions of velocity; and P 
and P a are even polynomials such that P( Q )(0) ^ and 
Pf a \(2) = 0. In the derivation, a scaling of Hp Hp in R 5 
appears as expected but it is eventually reduced to i? 4 by 
the incompressibility condition on /. Also noticeable are 
the vanishing contributions from pressure fluctuations in 
the torque du e to the symmetry of Q. 

In general, HpHp should not be invariant at large 
scales when R — > oo, as TpTp also appears to scale as 
R i . However, as with Loitsyanskii's integral, the behav- 
iors of f(s) and f a (s) when s — * oo must be taken into 
account in analyzing the scalings with R in Eq. (4) , since 
they can produce dive rgent in tegrals: for instance, with 
Saffman's f(s) ~ s~ 3 , Hp Hp ~ R 5 as first expected in- 



tuitively. For a generic f(s) 
it is thus found that 



behavior at s — > oo, 



4an. = lim 

it — »OG 



R 



m—4 ■ 



2R 



k I P(s/R)f(s,t)s 3 ds 
o 
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(5a) 
(5b) 



is a finite quantity for m < 4. No w, beca use in general 
fa(s) ~ / 2 (s) ~ s~ 2m at s — * oo, TpTp always scales 
as R for 1 < m < 4, and HpHp does become asymp- 
totically invariant at R — > oo for such values of m — or 
for < /x < 3 in spectral space. This proves the so-far 
empirical conditions for the permanence of large eddies 
reviewed in the introductory paragraphs. 
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At this point the significance of the invariant as given 
in Eqs. (5) requires some elaboration: (i) Although the 
physical invariant is always the same — the angular mo- 
mentum variance, — the expression of the invariant also 
depends on an external parameter m which is not spec- 
ified a priori: the long-range behavior of /(s) — or the 
infrared behavior of E(n). (ii) A given initial long-range 
behavior of /(s) will be preserved at later times accord- 
ing to Eq. (5a), i.e. m is invariant, if m < 4 — or /i 
if (j, < 3. (iii) Furthermore, if the evolution (by the 
Navier-Stokes equation) preserves the self-similarity of 
f(s,t), then kX m is also invariant according to Eq. (5b). 
(iv) As a corollary, because evolution does not affect the 
infrared spectrum below the integral length scale, and en- 
sures a quasi-equilibrium self-similar profile in the iner- 
tial range, it is necessary and sufficient to have initially a 
self-similar profile at large scales, such as /(s) ~ (s/A)~ m 
with m < 4 to ensure a kX m invariant, (v) For initial con- 
ditions at m > 4 the inert ial range dominates the integral 
in Eq. (5a), and H^H S J is not asymptotically invariant 
according to Eq. (4a). A universal behavior, i.e. inde- 
pendent of the initial infrared spectrum profile but not 
necessarily self-similar, must then appear and eventually 
converge towards m f=a 4 and n sa 4/3. This is the maxi- 
mum decay rate of turbulence. 

It must be stressed here that a rescaling of Loitsyan- 
skii's integral with the proper R power as for can- 
not produce the same invariants of HIT evolution be- 
cause different integration weights are involved — s 4 ds 
in Eq. (3) instead of s 3 ds in Eq. (5a). In particu- 
lar, the maximum decay rate at n w 4/3 cannot be 
predicted. This reflects a profound difference between 
the approaches of Loitsyanskii and Landau — despite the 
latter's incidental claim of equivalence 10 — which can be 
traced to alternative integrations of the stochastic equa- 
tion d t H^(R) = T^{R) (see Supplementary Methods, 
§1): Loitsyanskii's approach leads to d t H z H z — T Z H Z 
instead of Eq. (4a), and thus produces the two-point 
inpZe-correlation of the Karman-Howarth equation in- 
stead of the two-point quadruple-correlations of Eq. (4c) . 

Impulsive approach to Landau's i nvariant s. In the 
preliminary scaling analysis leading to H S2 H S2 ~ R 5 , the 
velocity field was implicitly assumed to display the struc- 
ture represented in Fig. 2a. Landau's result iLan. ~ R 3 
would actually correspond to Fig. 2b, where large eddies 
have no translational motion, and just rotational motion, 
each contributing by kX 2 instead of kR 2 to the angular 
momentum variance. 

The representations in Fig. 2, however simplistic they 
might appear, are actually well suited to carry out a scal- 
ing analysis and to derive the R dependence of the in- 
variants. Indeed the scaling behavior is marginally sen- 
sitive to the details of (i) the motions within the large 
eddies, contributing to lower powers of R, and (ii) the 
interactions and correlations between neighboring struc- 
tures. The last point stems from Saffman's projection 




Figure 2: Schematic representation of large turbulent eddies 
in HIT (shaded circles) moving randomly with typical veloc- 
ities ~ yk in two basic modes: a) translations or Saffman 
correlation, H 2 ~ R'kX 3 , b) rotations or Batchelor correla- 
tion, TP ~ R 3 kX s . 

procedure 6 which generates a kinematically consistent 
velocity field u from any given arbitrary impulse field 
i 

u = i — Vg with Aq = V- i, thus V- u = 0. (6) 

Accordingly, Fig. 2 represents impulse fields from which 
kinematically consistent velocities can be produced. 
Now, a crucial property of Saffman's procedure is that 
it strictly preserves the angular momentum of any spher- 
ical volume since, here again, the symmetry cancels 
the torque due to the "impulse pressure forces" Vg. 

The torque T/ 2 is not preserved by Saffman's proce- 
dure, but the i? 4 behavior of TpT^ is still recovered 
from a scaling analysis on the impulsive representation, 
as given in Supplementary Methods, §4. In particular, 
the differences of short range velocity correlations be- 
tween Figs. 2a and 2b do not impact the R 4 factor: in 
T/ 2 's expression, as the velocity field of the turbulent 
structures is just sampled at the surface df2, possible de- 
tails of the velocity correlation functions are lost. The 
basic R A scaling derived from Fig. 2a is then valid for all 
cases: according to Eq. (4c), only with strongly long- 
range-correlated velocities can the torque's scaling be 
changed — i.e. if f a (s) decays slower than s~ 2 or f(s) 
slower than s _1 . Rigorous and physically consistent con- 
clusions on both the scaling and the time dependence of 
Landau's invariants can thus be drawn from the present 
simple impulsive representation. 

As previously shown 6 , Saffman's projection gener- 
ally preserves the profile of the infrared spectrum: 
in the canonical case of Fig. 2a first considered by 
Saffman 6 , random translations produce E(k) ~ n 2 . At 
the same time however, it can dramatically modify the 
long range profiles of the two-point self-correlations / 
of i and it 6 : again in Fig. 2a, an arbitrary i represent- 
ing random long-range-uncorrelated translations — thus 
with transcendentally vanishing fi(s) — eventually yields 
a long-range-correlated it with an algebraically vanishing 
f(s) ~ s~ 3 — reflecting the laminar and dipolar charac- 
ter of long range correlations steming from Eq. (6) . This 
last result substituted in Eq. (4b) recovers the proper 
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Figure 3: Impulse representation of large eddies undergoing 
so-called "backscattering" in Batchelor type HIT. The process, 
albeit slower, also affects "big frozen eddies", much larger than 
the integral length scale. 

behavior I^ an , ~ R 5 as estimated from the impulse field 
i. 

The E(k) ~ k 4 infrared profile first considered by 
Batchelor & Proudman corresponds to a purely rota- 
tional impulse field as in Fig. 2b: it is indeed found 
Ihan. ~ R 3 in this case, both with the impulsive scaling 
analysis and with Eq. (5a). Since m > 4, the expected 
non-self-similar evolution is obtained where, akin to the 
collisionnal processes at molecular scales, the turbulent 
energy is in part redistributed on translational degrees 
of freedom with lower m < 4 — a well-known "backscat- 
tering" phenomenon illustrated with an impulse image 
in Fig. 3. H owever, the typical turnover rate at scale R 
is given by (Hp Hp) 1 / 2 / '(pR 5 ) ~ R~ 2 and at very large 
scales the redistribution is practically frozen. 

Saffman's impulse field appears as most likely when 
considering turbulence generated by experimental set ups 
(as grids) devoid of specific devices which could constrain 
correlations such as Batchelor's 6 . However, (i) in many 
numerical simulations, initial conditions are implicitly of 
Batchelor type, and (ii) any experimental or numerical 
setup has an upper bound for R and produces Batchelor 
fields at these largest scales. Other types of initial con- 
ditions seem unrealistic or highly specific: for instance 
E(k) ~ k for k — ► ensures full self-similarity of the 
whole spectrum down to dissipative scales 11 , but no gen- 
eral physical mechanism is known to produce such a con- 
straint. 

Invariants in layer and tube. The impulsive approach 
to turbulence invariants has the important property of 
being easily extended to the less trivial situations of self- 
similar layers and tubes of turbulence. Only the z com- 
ponent of angular momentum along the symmetry 
axis is now preserved: the previous derivation for HIT 
is simply adapted to these cases taking the integration 
volumes sketched in Fig. 4. Details are provided in the 
Supplementary Methods, §4, but the basic steps and as- 
sumptions are quickly summarized below. 

Because of self-similarity, typical layer thickness or 
tube diameter scale as A. As in HIT, for each motion 
mode (Saffman, Batchelor, etc.) of each geometry, fac- 
tors in Hp Hp come from velocity fluctuations, typical 
distance to axis, and integration domain over correlated 
and uncorrelated distances — below and beyond A. With 
vanishingly small contributions of velocity fluctuations 
in the surrounding inviscid laminar fluid, the torque is 
dominated by the parts of f2's surface which intersect the 




Figure 4: Volumes Q(R) (framing lines) displaying asymptot- 
ically invariant variance of angular momentum for turbulent 
layer, tube, and spot (shaded volumes). 

turbulen t zone. Then, comparing behaviors as R — ► oo, 
Hp Hp is also found to be asymptotically invariant for m 
below some maximum value. Finally, Saffman's projec- 
tion procedure applies as well: impulse is confined to the 
turbulent zone, velocity extends into laminar regions, but 
angular momentum is preserved by the projection. Re- 
sulting invariants and evolution exponents are collected 
in Table I for Saffman impulse fields, and for correlations 
corresponding to maximum decay rates. The strength 
of the impulsive approach is now apparent as all these 
new non-trivial results are readily obtained "on the back 
of the envelope," without any detailed knowledge of the 
complex correlation functions that describe the velocity 
field in these anisotropic and inhomogeneous flows. 

As reported in Table I, Batchelor and Saffman invari- 
ants happen to coincide for the tube of turbulence be- 
cause both a re produ ced by identical fcA 2 as point con- 
tributions to Hp Hp. The case of a spot of turbulence 
is also considered but, in contrast to the three previ- 
ous situations, the angular momentum cannot fluctuate: 
mean angular momentum due to the residual overall spin- 
ning then provides the invariant in Table I. The turbu- 
lent tube is also compatible with mean overall spinning — 
centrifugal forces remain finite — with the corresponding 
invariant in Table I. The relaxation of turbulence for spin- 
ning zones is constrained by these invariants for as long as 
mean angular momentum dominates other terms — such 
as angular momentum variance in the tube case. 

Besides HIT, experimental and numerical results 
appear to be available for the turbulent layer 
only, a situation also known as the incompressible 
Richtmyer-Meshkov instability (IRM) at vanishing At- 
wood number 12 ' 13 . The value of the IRM self-similar 
growth rate is still a matter of debate due to many per- 
turbation sources affecting accuracy (limited time span 
and space extension, poorly defined or controlled ini- 
tial conditions, etc.). However, up to reported uncer- 
tainties, values measured for 1 — n/2 do fall within the 
expected [2/7,1/3] range, with results from experiments 
and from numerical simulations closer respectively to 
the upper 1/3 and lower 2/7 bounds. This is consis- 
tent with the corresponding initial conditions: expected 
Saffman type in experiments and explicit Batchelor type 
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Table I: Turbulent relaxation in geometries of Fig. 4: invariants and evolution exponents (S. & B. = Saffman & Batchelor). 



Geometry 
Invariant type 




HIT 


Layer (IRM) 


Tube 




Spot 


Saffman 


Maximum 


Saffman 


Maximum 


S. & B. 


Spin 


Spin 


m (invariant k\ m ) 


3 


4 


4 


5 


7 


6 


8 


n (decay exponent of k) 


6/5 


4/3 


4/3 


10/7 


14/9 


3/2 


8/5 


1 — n/2 (growth exponent of A) 


2/5 


1/3 


1/3 


2/7 


2/9 


1/4 


1/5 


C £ 2 (e model constant) 


11/6 


7/4 


2 


15/8 


19/10 


2 


2 



in simulations — in the latter, calculations start from su- 
perpositions of random interfacial perturbations in the 
laminar linear regime, hence with vanishing . A pre- 
vious theoretical prediction based on an implicit assump- 
tion of Saffman type turbulence also yielded 1 — n/2 = 
1/3 as expected, but with a completely different approach 
which could be extented in 2D 14 . 

The present predictions for different flow topologies — 
of higher practical relevance than HIT — provide guidance 
in the modeling of the equation of turbulent dissipation 
e 15 . For instance, e is usually modeled as in the k- 
e model which, for the four situations considered here, 
can be approximately but accurately reduced to the bulk 
ODEs 16 

d t K=-d(d t A/A)K-E (7a) 
d t E = -d (d t A/A)E - C E2 E 2 /K (7b) 

where K, E, and A are the bulk averages of k, e and 
A over the turbulent zones, dtA/A terms describe the 
dilution effects in place of the flux terms of the PDEs, and 
d = 0, 1, 2, or 3 for HIT, layer, tube or spot respectively. 
Self-similar solutions of Eq. (7) match the invariant based 
solutions if the model constant is 

C e2 = 3/2 + (1 + d/2)/(m - d) (8) 

yielding the values in Table I. 

Inspection of Table I shows that a single value of 
C £ 2 — usually set to 1.92 15 — can never capture the three 
Saffman cases at once (let alone the others), possibly ex- 
plaining some practical difficulties such as up to 15% er- 
rors in growth rates when trying to simultaneously cap- 
ture plane and cylindrical jets with a single model 15 . Sig- 
nificant improvements could thus be contemplated when 
adding model equations aimed at "sensing" the dimen- 
sionality of the turbulent zone d, and at "defining" the 
behavior of angular momentum variance at large scales 
to. Anyhow, these correction will remain around 15% 
since the influences of m and d in Eq. 8 almost compen- 
sate each other when varying the flow topology. Inciden- 
taly, this noticeable property of Eq. 8 explains in part 
the often disregarded fact that most models are robust 
with respect to the effective Knudsen number of turbu- 
lence: models can and are applied to situations where A 
is similar to the characteristic flow size L without major 



malfunctions, although they are calibrated through G e2 
to capture HIT where A -C L. 

The generally accepted perception of the significance 
and status of e in turbulence modeling 10 should be re- 
assessed in view of the present results. The initial defini- 
tion as the mean small-scale viscous dissipation VUiJuJJ 
is now often abandoned in favor of the spectral flux, 
which assumes some form of spectral quasi-equilibrium 
whereby E(k) oc £ 2 / 3 k -5 / 3 in the inertial range. As 
shown here however, turbulence decay is generally con- 
trolled by the infrared spectrum regardless of the detailed 
profile in the inertial range: the very same results would 
hold even with non-Kolmogorov spectra — which do ap- 
pear for instance in the presence of helicity, buoyancy or 
electromagnetic effects. In fact, e is solely defined as the 
time derivative of k, the energy in the energy containing 
range. Thanks to the assumption of global self-similarity 
of the energy containing range (infrared + inertial), an 
integral length scale A can be defined, and then e is con- 
strained to match fc 3 / 2 /A. 

It should not be considered as surprising that the self- 
similar decay of turbulence is totally independent of the 
detailed processes within the inertial range. Indeed, this 
situation is comparable to the perfectly inelastic colli- 
sion of two bodies: the overall conservation of momen- 
tum fully defines the final state, regardless of whatsoever 
details on the mechanisms which have actually removed 
or dissipated the excess energy. 

Conclusions. The present work has provided a rigor- 
ous answer to the long-standing issue of the permanence 
of large eddies and their relationship to turbulence decay, 
using angular momentum invariance at large scales. This 
approach has the advantage of being both extendable to 
less trivial cases of turbulent layer, tube and spot, and 
reducible to a calculationally simple impulsive picture. 
Beyond their intrinsic value and novelty in the theoreti- 
cal study of turbulence, the resulting predictions on decay 
rates open up new approaches and understanding for im- 
proving the modeling of turbulent dissipation in applied 
simulations. 
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SUPPLEMENTARY METHODS 
1. Landau's angular momentum invariance at large scales 

One starts with the equation of fluctuating velocity in a zero mean velocity field (uniform constant density and 
vanishing viscosity are assumed) 

pd t Ui + p(muj)j + pj = 0. (9) 
Angular momentum fluctuations along z integrated over a fixed volume Q then yield 

/ CzijU [pdtUj + (pujU k ). k +p,j] d 3 r = 0, (10) 
Jn 

or after elementary algebra 

d t Ff=Tf, (11) 
where the fluctuating angular momentum and torque are 

Hp (f) = f , .,,r,p„,A\ '-r. (12a) 
Jn 

T ^(t) = - £zijri(puju k + Sjkp) <7fcd 2 r, (12b) 
Jan 

<x being the unit vector normal to the surface, with components <7j. For integration volumes f2 considered in this 
work, the orientation of the normal vector element <xd 2 r cancels the z-component of the torque produced by pressure 
p. Pressure is thus disregarded in all the following. 

A time evolution equation for the variance of Hp can then be obtained by a direct substitution of the torque 
expression, and accordingly 



d t Hp(t)Hp(t) = 2 Tp(t)Hp{t) = -2 // WVi'i-ri A>KW^(r> t d 2 rdV. (13) 



anxn 

In this way, the two-point, triple space-correlation functions of velocity can be introduced, in close similarity with the 
usual Karman-Howarth equation in volume integrated form. Now, this derivation, although rigorous and exact, is not 
in the spirit of Landau's approach, which hinges critically on the argument that the torque becomes asymptotically 
negligible at large scales, and the angular momentum asymptotically constant (one is a surface integral and the other 
a volume integral) . In this case one should instead derive an evolution equation of the angular momentum variance 
by integrating the stochastic differential equation (11) in the weak collision (or friction) limit — just as for Brownian 
motion. For asymptotically large fl the correlation time of Hp(t) is much longer than that of T' z ?(t), and one can 
thus write 



d t Hp(t)Hp(t) = 2 Tp{t) / Tp{t')dt> « 2 / T»(t)Tf(t - t')dt' = 2 T?(t)T?(t) r(t), (14) 

Jo Jo 

where the torque correlation time is defined and approximated as 

i r°° 

T{t) = Wmml vwn- (15) 

Here, it was assumed that the spinning associated with the angular momentum is small enough so as to produce 
negligible shear at the boundary: this avoids far-from-average velocity fluctuations which could significantly change 
the torque estimate and even p roduce a correlation between tor que and angular momentum. 

In general, expanding TpTf shows that the time derivative of Hp Hp now depends on two-point, quadruple space- 
correlation functions of velocity. Landau's approach is thus physically distinct from Loitsyanskii's, although eventual 
invariants ma y — and shou ld — coi ncide in regimes where their respective approximations are compatible. Explicit 
expressions of Hp Hp and TpTp as integrals of velocity space-correlation functions are given in the next two sections 
for the simple case of HIT. 



8 



2. Angular momentum variance of a sphere in HIT 

The angular momentum of a given spherical volume Q of radius R around origin in a HIT field is given by Eq. (12a). 
Although it cancels on average, its average square norm does not and is given by 



jktij'k'rjr'j, p z u k {r)u k ,{r') d 3 r dV 



P 2 JJ ^Sjj'Skk' - Sjk'Sj'k) rjr'j, u k (0)u k ,{r> - r) d 3 r dV 
\ p 2 k [[ 2 (r jr >S kk , - r k ,r' k ) (g{s)5 kk , + [f(s) - g(s)]^) 



d 3 r d 3 r' 



n 2 

= \p 2 k 

JJn 2 

= \p 2 k 
6 JJn 

= \p 2 k 
A JJn 2 



J( s ) +g ( s )]r k r' k -[f(s)-g(s)} 



[f(s)+g(s)} [x 2 -s 2 /4]-[/(s)- 5 (s)] 



/(*) 



(x ■ sf 



g(s) 



n (X • S) 2 o . 

x + * s-i- - s 2 /2 



-« a /4 
d 3 sd 3 j;, 



d 3 sd 3 a; 



(16) 



where homogeneity and isotropy have been taken into account, with s = r' — r, x = (r' + f)/2, and / and g being 
the longitudinal and transverse correlation functions. The integration domain J? 2 is represented in x and s variables 
as x 2 ± x ■ s + ,s 2 /4 < R 2 , so using polar coordinates of x with respect to s (with c~x- s/xs), it is found 



d 3 sd 3 x = 8tt 2 



n- 



2R ( j-R-s/2 pi 
J-l 



x dc dx 



v /R 2 -s 2 /4 p(R 2 ~x 2 -s 2 /i)/xs 



R-s/2 



-(R 2 -x 2 -s 2 /4)/xs 



x 2 dcdx s 2 ds, (17) 



and thus, after lengthy but straightforward calculations 



Hp Hp = R 5 —p 2 k 
45 



f 2R / 9s 3s 2 \ / s \ 3 



16- 



13s 21s 2 21s 3 



2R 2R 2 8R 3 



s \' 
1 -2r) 



g(s) 



s 2 ds. (18) 



Using the incompressibility condition g(s) = /(s) + s/'(s)/2 to eliminate g(s) and performing integrations by parts 
on s™/'(s) terms to eliminate f'{s), eventually yields 



Hp Hp = R 4 —p 2 k 
3 



2R 



2R 2 J V 4i? 2 J 



f{s)s 3 ds. 



(19) 



It is to be noticed that the R 5 J fs 2 ds and f(2R) terms (the latter from partial integrations) both cancel exactly — not 
just as R — > oo, — leaving i? 4 as leading term whenever J fs 3 ds converges at oo. 



3. Torque variance on a sphere in HIT 



The torque on a given spherical volume Q of radius R around origin in a HIT field is given by Eq. (12b). Although 
it cancels on average, its average square norm does not and is given by 

WW= I! eijkeij'vrM, p 2 u k {r)ui{r)u k > {r')u v {r') aia[, d 2 rdV 
JJ(dny 2 

= R 2 P 2 (( (aja'jSkv - a y a' k ) u k {0)ui{Q)uy \r' - r)u v {r' - r) cna[, d 2 r dV, (20) 
J J (an) 2 

where now <r = r/r = r/R is the unit vector normal to the sphere surface dfi. Here as for Hp Hp in the previous 
section, it is necessary to provide the general two-point quadruple velocity correlation tensor, which in HIT can be 
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decomposed as 

2 



— ) u fc (0)wi(0)u fc /(s)w;/(s) = fi(s) 



SkSlS k >Sli 



Sk'Sl' 



SkSl 



+ .91 0) S k i — n- + —j- x , 



SfcSiSfc'S/' 



A-'/' 



SfcSfc' SfcS;' 

S— 0«' H 7T- 



92 (s) 



+ h 2 (s) ( Skk'Sw + 5 k i> x ik> 





SlS k ' 


X j_ SlSV 


, , SkSiS k >Si> \ 


s 2 


OkV H 5- 


■**» 4 S 4 J 


SkSl 


Sk'i' - 




! SkSiSk>sv \ 


s 2 


Ofc; 9 - 


s 2 ) 


SkSk 

s 2 


- 5 W - 




sis k > , stsy 

9 <Jfci' 9 <Jfcfc 



SkSl 



k'V 



2— + 25ki k J 
s s 



(21) 



corresponding to the five different symmetrized and orthogonalized correlations tensors along zz\zz, zz\xx 7 zx\zx, 
xx\yy 7 and xy\xy when s is along z. The correlation functions have the following elementary properties 



32(0)=. 9i (0), 
fta(0) = hi(p), 
/i(s), 5 i(s),/i 1 (s) > 0. 



lim f 1 (s),g 1 (s),h 1 (s) = 1, 

s — >oo 

lim g 2 (s),h2{s) = 0, 



Taking into account that <x and it' are unit vectors, that s = r' — r = R{&' — &), and defining <r • a' 
1 — s 2 /(2i? 2 ), the various contractions of the elementary tensors yield 

, I s I x SkSlSk'SV i c 2 - 1 
((TjCTjdkk' -Ok<V k ) "4 a \°V = A ' 



f * k >\ fx Sk ' Sl ' i SkSl x \ i n 

{0-j0 3 bkk' - Vk'Vk) y °kl 1" — J" <W J (7;0> = 0, 



{(Jj(j'j5kk' - crk'cr'k) 



SkSk' 



, S k Sl' SjSk' , 

I' H 9— <W H 5— dfep 

s z 

((TjO-jSfck' - Vk'&k) 8kl X k'l'(7l<Tl' 



5c 2 - c -2 



Therefore 



a 



R- 



(dO) 2 



fi(s)-2g 1 (s) + h 1 (s)) + 



3c 2 - 5c 



92 (s) 



0, 
3c 2 

I- 5c 



(22a) 
(22b) 
(22c) 
cos6> = 

(23a) 
(23b) 

(23c) 

(23d) 
(23e) 



d 2 rdV. (24) 



Using polar coordinates relative to r for instance, so r ■ r' = r 2 cosO, and substituting cos# = 1 — s 2 /(2R 2 ) 



TpTP = 2R' 



s 

2R 2 J 
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i? 2 



h 9 (s) 



s 2 \ 



s 

R 2 



1 



R 2 J 

— ) 

AR 2 J 



92(s) 



)[fi{s)-2g 1 (s) + h 1 (s)]\sds. (25) 



Here again the incompressibility conditions could produce some simplifications. It is necessary however, to first 
obtain these conditions on the two-point quadruple correlation tensor in Eq. (21). The tensor is purely rotational 
in any of its components, and by symmetry, a zero divergence with respect to the first component suffices to ensure 
incompressibility. Now, from its general form in Eq. (21), substituted by monomial tensors T£ lk , v and functions f a 
as combinations of /1 to h 2 , and using s ^ = Sk/s, the incompressibility condition can be written 



E f«( s ) T Mk>v(s) B lwv = E S» t T m>i'Bwv + f a {s)T% im , tk B im , = 0, 

> K a 



(26) 



where Bik'v is any arbitrary tensor. The decomposition of the third-order tensors Bik'i> in irreducible representations 
taking into account the symmetry properties shows that three scalar conditions are eventually obtained, corresponding 
to the three following orthogonal tensors Bik'v 



sisk'sv 



81 x 

— Ok' l' 

s 



SlSk' Si' 



x Sk ' 



SlSk'Sf 



(27) 
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or equivalently, to the three independent tensors 

SlSk'Sl' 



££ 
S 



Sk' 

s 



The table of the different contractions (s&/ ' s)T klk , ir Bik>i' and T klk , v k Bik'v is thus obtained: 



1 klk'l' 








s 3 


Sl' 


s 


<w 


Sw 


Sfc' 

s 


SkSlSk'S V 

« 4 
a 




SlSk'Sl' 




1 


2 

s 


1 


2 

■S 


1 


2 

■S 


, Sk'Sl' 
5 kl 2 

SkSl 

s z 


S 


2 9 Ofc'i' 

s^ 


SlSk'Sl' 



1 
1 

1 



2 

s 


1 






6 

s 


1 
1 


2 

■S 

2 

s 


SkSk' r- 

s 1 

Sk&l' c 

s^ 




2 -4- 

s^ 

S// 

2 -9- <W 

s^ 




1 
1 


2 

s 
2 

s 


1 
1 


2 

s 
2 

s 


3 
1 


(i 

u 

S 

2 

s 


9 OH' 

s z 


Sk' 

—z-011' - 
s l 


- -7 Ok'l' - 

s z 


SlSk'Sl' 
S 4 


1 





1 


2 

s 


1 


2 

s 


SlSl' 

9 Ofcfc' 

s z 


s l' p 

-o-Oife' " 

s z 


-%8 k ' V - 
s z 


SlSk'Sf 
S 4 


1 





1 


2 

s 


1 





Sklh'l' 









1 





3 





1 





0~kk'0~ll' 









1 





1 





3 





hi'hn 









1 





1 





1 






Hence the three conditions 



sf[/2 + f x - gi - 2. 92 = 0, 
sg[/2 + g 1 +g 2 -h 1 +h 2 = 0, 
sg'2 + 3.g2 + gi - hi - h 2 = 0. 



(28) 



(29a) 
(29b) 
(29c) 



To reduce the R 4 scaling of TpT t n in Eq. (25), the leading integrand term proportional to (3h 2 — g 2 )sds should be 
canceled by combining the incompressibility constraints and performing integrations by parts. Taking coefficients a, 
/?, and 7 for the three constraints, this leading term is written as: 

[(3/i2 - 92) + a{sf[/2 + /1 - gi - 2g 2 ) + P{sg'j2 + gi + g 2 -hi + h 2 ) + 1 (sg' 2 + 3g 2 + gi - hi - h 2 )\ sds 

= [ - (a - 7)01 - (2a - p - 7 + 1)<? 2 - (/? + 7)^1 + {ft - 7 + 3)^] sds, +d[s 2 (af 1 + 0g x + jg 2 )] . (30) 

In general, no combinations of a, (3, and 7 can be found to simultaneously cancel the four integrands on the right 
hand side. The scaling in R 4 as given in Eq. (25) is thus completely general, and the incompressibility conditions just 
provide the possibility of expressing the leading term as a function of only one of the four correlation functions <?i , 
<?2, hi, or h 2 — a substitution of minor interest. 



4. Scalings of HpHp and TpTp in the basic C x configurations 



1. Generic scaling approach. As mentioned in the 
main text, the impulsive picture of the velocity fluctu- 
ations can accurately provide scaling laws for H^ 2 H[ 2 
and TpTp in HIT, and in layer and tube of turbulence. 



This is ensured by the basic property of the Saffman pro- 
cedure which produces kinematically consistent velocity 
fields out from arbitrary impulse fields, while preserving 
at the same time the angular momentum . 

The case of HIT is represented in Fig. 1. Within struc- 
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tures of typical size A, the velocity field is modeled as fully 
correlated. In the Saffman case, this produces c ontribu- 
tions from each structure to the double integral in Hp Hp 
scaling as ~ kR 2 X 6 . The scaling is reduced to ~ kX 2 X 6 in 
the Batchelor case where motions are purely rotational, 
since opposite contributions from within a given struc- 
ture cancel in first order of R. Between structures, the 
velocity field is now modele d as fully uncorrelated: the 
double sum on fi in Hp Hp multiplies the structure con- 
tribution by ~ (R/X) 3 , the typical number of structures 
in volume fi. 

This gene ral app roach to estimating scalings is simply 
adapted for TpTp and for the layer and tube of turbu- 
lence. The various intermediate contributions and final 
results are collected in Table II, but some peculiarities 
require further comments as given below. 

2. Relationship with rigorous calculations for HIT. 
As already menti oned in the main text, the rigorous cal- 
culation in §2 of (iff) 2 for HIT gives Eq. (19), with a 
leading contribution in R 4 under incompressibility con- 
straints which are not fulfilled by the impulse fields con- 
sidered here. A R 5 behavior is actually found in the 
general case in Eq. (18), compatible with the Saffman 
case. A R 3 behavior could be recovered with a specific 
oscillating profile of the correlation functions / and g, as 
would be expected in the Batchelor case. 

3. Scaling of TpTp. The e stimatio n of TpTp follows 
the same general lines as for Hp Hp. For the Saffman 
case of HIT in Fig. la, contrib utions of each structure to 
the double surface integral in TpTp scale as ~ k 2 R 2 X i . 
Now, at variance with the estimation of Hp Hp given 
above, this scaling is conserved in the Batchelor case of 
Fig. lb and must not be reduced to ~ k 2 X 2 X 4 : opposite 
contributions from within a given structure are not inte- 
grated over its full volume but on a randomly oriented 
slice only, and therefore the cancellation of the torque in 
first order of R is incomplete. Finally, in all c ases of im- 
pulse fields, the double sum on dfi in TpTp multiplies 
the structure contribution by ~ (R/X) 2 , the typical num- 
ber of structures intersecting surface dfi. Therefore, the 
TpTp scaling is independent of the type of the correla- 
tions in the impulse field, and sets the maximum scaling 
exponent of Hp Hp compatible with asymptotic invari- 
ance at R — > oo. As already mentioned in the main text, 
although the Batchelor impulse field does not produce 
an invariant, it should still be considered here for its rel- 
evance in experiments and simulations, where physical 
boundaries to the flows are always present at som e large 
scale. This important scaling property of TpTp is also 
found in the layer and tube of turbulence, and Table II 
thus lists the corresponding invariants of maximum ex- 
ponents. 

4- Basic scalings in configurations other than HIT. 
The major difference between the different flows con- 
sidered here comes from the double integral over fi in 



HP Hp: proportional to (R/X) 3 , (R/X) 2 , and R/X for 
HIT, layer, and tube respectively. An interesting fea- 
ture is also to be noticed on the point contribution to 
Hp Hp for the tube of turbulence: the typical distance 
to the z axis of any given point in the turbulent zone is 
A, instead of R for the HIT and layer of turbulence. Ac- 
cordingly, all the possible types of impulse field s — here 
the Saffman and Batchelor cases — yield identical Hp Hp 
scalings. The evolution of the tube of turbulence is thus 
more universal as it is independent of the long-range cor- 
relations in the initial conditions. 

5. Spinning tube of turbulence. The integration vol- 
umes fi used for deriving asymptotic invariants are all of 
Coo symmetry and are thus also invariant by rotations. 
It is therefore natural to consider non-vanishing mean 
velocity fields producing global spinning of the turbulent 
zones. However, aside from the spot of turbulence exam- 
ined below, the tube of turbulence is the only flow config- 
uration where a homogeneous spinning will not produce 
divergent centrifugal forces. A spinning tube of turbu- 
lence could see its turbulence relaxation contro lled by 
the conservation of the angular momentum itself, Hp in- 
stead of Hp Hp . This somewhat different situation im- 
plies that (i) the mean velocity field is non zero in the 
turbulent zone, (ii) the flow evolution produces an equili- 
bration between the mean kinetic energy contained in the 
angular momentum and the turbulent kinetic energy — a 
reasonable assumption, — and (iii) the mean angular mo- 
mentum dominates its variance, otherwise the evolution 
would b e constrained by the large scale invariance of 
Hp Hp. The scaling in Table II is then found, given 
as a function of Hp x Hp for consistency with the other 
results. 

6. Spinning spot of turbulence. Because the extension 
of a spot of turbulence is finite in all directions, it is not 
possible to predict its evolution from angular momentum 
variance at large scales. More directly, the angular mo- 
mentum itself is invariant for the volume fi in Fig. 3. As 
for the spinning tube of turbulence, this implies (i) a non- 
zero average velocity field, (ii) an equilibration between 
the mean and turbulent kinetic energies, and (iii) a non- 
negligible angular momentum, otherwise the evolution of 
the corresponding quadrupolar — or higher multipolar — 
spot would not be constrained by any obvious invariant. 
T he sc ali ng in Table II is then found, again as a function 
of Hp x Hp. 

7. Axes other than z for angular momentum conser- 
vation. In principle, for the geometries other than HIT 
of Fig. 3, one c ould con sider the asymptotic invariance 
of Hp Hp — not Hp Hp — for a spherical volume fi, in- 
stead of slab and rod shaped volumes. Pressure forces 
would still not contribute to the torque, but the fringe 
field of velocity fluctuations in the laminar regions would 
contribute substantially and in a complex way to both 
the angular momentum and torque, as the laminar re- 
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Figure 5: Schematic representation of the build up of velocity 
fluctuations and correlations in the laminar surroundings of a 
turbulent zone (layer or tube). Arrows represent the potential 
field contributions of turbulent structures affecting a given 
point of the laminar zone. Squares represent typical domains 
of correlated fluctuations. 

gions would represent the largest volume fractions of the 
sphere. The corresponding invariance conditions could 
then actually not constrain the evolution of the turbu- 
lence zone, although they would b e compatible with the 
present findings based on Hp Hp. 

8. Contributions to Hp Hp from laminar zones of fl. 
As shown in Fig. 3, the domains fl considered for angular 
momentum conservation at large scales extend somewhat 
beyond the turbulent zones. In these laminar regions, ve- 
locity fluctuations exist due to the long range influence of 
the motions in the turbulent zone, although they decay 
quickly with distance. These velocity fluctuations could 
contribute to Hp ', despite the fact that Hp is preserved 
by Saffman's projection procedure in Eq. (6): indeed, 
Saffman's projection only deals with a pressure gradient 
correction to the impulse field, but does not involve trans- 
port which is the source of torque in Eq. (12b). Now, as 
velocity in the laminar region can be described by a po- 
tential flow whereby u = V0, the contribution to Hp 
of a given volume fli in the laminar zone of fl with the 
same C x symmetry is 

H "* = / £*ij r ipUj d 3 r 




= / e z ijnp<f>(Tj d 2 r = 0. (31) 
Jdn t 

The final cancellation comes from the symmetry of fli, 
sorXff has no component along z (<x is the unit vector 
normal to the surface dflg). It also assumes that flu is 
simply connected, or is just connected but does not carry 
any net circulation of u around singularities (as could be 
the case for the tube of turbulence) . 

9. Contributions to TpTp from laminar zones of dfl. 
As already mentioned, angular momentum cancels in the 



laminar regions but can be transported across these re- 
gions to and from the turbulent zone. It is thus necessary 
to evaluate the magnitude of the torque contributions on 
the side surfaces of the domains fl represented in Fig. 3. 
As with the contribution from the t urbulen t zone esti- 
mated above, the double integral in TpTp will be de- 
fined by the typical velocity fluctuations at each point, 
summed over the surface in correlated or uncorrelated 
ways. 

Let d be the distance of the surface of fl to the tur- 
bulent zone. The velocity fluctuations at distance d are 
given by the potential field produced by the turbulent 
zone in the laminar region, which scales as (X/d) 3 Vk or 
(X/d) A Vk for a single structure of dipolar or quadrupo- 
lar character respectively (Saffman or Batchelor cases). 
As sketched in Fig. 5, the effective field at distance d is 
the superposition of uncorrelated contributions coming 
from structures within a typical radius d: this yields the 
various point contributions as listed in Table II. Since 
at distance d the field is produced by turbulent struc- 
tures within a typical radius d, velocity fluctuations are 
also correlated over a typical distance d: this yields the 
various surface integration factors as listed in Table II. 
Combining the different factors , the re lative magnitudes 
of the laminar contributions to TpTp are then found as 
listed in Table II. 

In order to preserve the sam e invaria nce c onditio ns as 
defined by the scalings in R of Hp Hp and TpTp , it is 
necessa ry to make the relative laminar contribution to 
TpTp decrease with R. Upon inspection of results in 
Table II, this is always possible in all cases by making d 
vary with R, for instance as d oc i?"^ 1- "-* with 1/5 < 
a < 1. This justifies, a posteriori, the requirement that 
fl extends somewhat into the laminar region. 

10. Contributions to TpTp from laminar zones of dfl 
in spinning tube and spot. The derivation of the previ- 
ous paragraph must be somewhat adapted to the cases 
of spinning tube and spot of turbulence. The fluctua- 
tions in the surrounding laminar fluid due to the rota- 
tion of the turbulent zone come from the irregular shape 
of the fluid domain entrained by the rotation. The po- 
tential flow around a spinning spot of turbulence would 
then produce velocity fluctuations decaying at least as 
(A/d) 4 , as for a Batchelor structure, faster decays being 
possible for shapes of rotating fluid domains of higher 
multipolar structure. Hence the results for the spot of 
turbulence listed in Table II. It should be noticed that, 
for the spot, the torque on the surface dfli in the lam- 
inar domain actually rep resen ts the full torque, and the 
asymptotic invariance of Hp x Hp is now obtained for 
d — ► oo. The results for the spinning tube of turbulence 
are directly adapted from the spot case, assuming un- 
correlated phases and shapes of the spinning sections of 
typical size A over length R. 



Tabic II: Magnitudes of leading terms in R to Hp Hp and TpTp for the various flow configurations and impulse field types considered here, and corresponding torque 
contributions from laminar zones (TZ = turbulent zone, LS = laminar surroundings). 



Geometry 


HIT 


Layer (IRM) 


Tube 




Spot 


Invariant type 


Saffman Batchelor 


Saffman Batchelor 


Saffman Batchelor 


Spin 


Spin 


Point contribution to HpHp in TZ 

Sum over volume in TZ 

Hp Hp scaling in TZ 

Point contribution to TpTP in TZ 

Sum over surface in TZ 

TpTP scaling in TZ 


kR 2 k\ 2 
xA 6 x (R/Xf 
R 5 kX 3 R 3 kX 5 
k 2 R 2 
xA 4 x (R/X) 2 
R 4 k 2 X 2 


kR 2 kX 2 
xA 6 x {R/X) 2 
R 4 kX 4 R 2 kX 6 
k 2 R 2 
xA 4 x (R/X) 
R s k 2 X 3 


kX 2 kX 2 

xA 6 x (R/X) 
RkX 7 RkX 7 
k 2 X 2 
xA 4 
k 2 X 6 


kX 2 

xA 6 x (R/X) 2 
R 2 kX (i 
k 2 X 2 
xX 4 
fc 2 A 6 


kX 2 
xA 6 
fcA 8 

(A/d) 10 fc 2 A 6 


Point contribution to TpTp in LS 
Sum over surface in LS 

TpTp: scaling of LS to TZ contributions ratio 




(X/dfk 2 R 2 (X/d) 12 k 2 R 2 
xd 4 x (R/d) 2 
X 5 R/d 6 X 9 R/d 10 


(X/d) u, k 2 d 2 (X/d) 14 k 2 d 2 
xd 4 x (R/d) 
X 4 R/d 5 X 8 R/d 9 


(X/d) 14 k 2 d 2 
xd 4 x (R/d) 
X s R/<? 


(X/d) 16 k 2 d 2 
xd 4 



